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Manifestations of the axial anomaly in finite temperature QCD* 
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We compute the flavor singlet meson correlators and screening masses in quenched and Nf = 2 QCD at 
Nt = 8. The consequences of our results for the realization of the Ua (1) symmetry at finite T are discussed and 
an interpretation of our measurements in terms of the behaviour of the low lying fermionic modes is proposed. 



Although the QCD chiral phase transition has 
been studied for many years, relatively little is 
known so far about the realization of the axial 
U(l) symmetry in finite temperature QCD. The 
limiting cases are understood however: at very 
high temperature the Ua{1) symmetry can be 
considered as "effectively restored" (since the ef- 
fects associated with the axial anomaly are small 
in a dilute gas of instantons) whereas at zero tem- 
perature a sizeable explicit breaking of the L0t(l) 
symmetry is essential to our understanding of me- 
son spectroscopy. The question therefore arises as 
to how the transition between these two regimes 
is realized and how this relates to the restora- 
tion of the SU(Nf) chiral symmetry. In order 
to investigate this problem, we have computed 
mesonic screening correlators at various tempera- 
tures above and below the phase transition. Con- 
sidering 2 flavors of valence quarks, we do this in 4 
channels corresponding to the 7?, a, 5 and rf . Chi- 
ral symmetry restoration would then induce de- 
generacies represented by the following grouping 
of states: (tt, a) and (6,r]'), while J7a(1) restora- 
tion would imply: {<r,rf) and (tt, S). 

Our lattice computations are done using stag- 
gered quarks. The operators representing the w, 
a, S and 7/ are respectively 75<8><!;5, I® I, I®^5 and 
75 ® J. The first two being local operators and the 
last two 4-link operators. In order to minimize 
the effect of flavor symmetry breaking (which can 
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be substantial at the values of (3 currently used) 
it is desirable to limit comparisons to one of the 
two categories. This strategy allows a study of 
the 7?, a, S system with minimal lattice artifacts 
if one replaces the 5 correlator by the connected 
part of the a correlator {a conn ). This being said, 
the computations involving the rf operator are 
also quite important in exposing the underlying 
physics, since (through the Atyah-Singer index 
theorem) 75 ® / serves as an indicator of topo- 
logical activity in the QCD vacuum. All our cor- 
relators (connected and disconnected) were com- 
puted using a U(l) noisy estimator following the 
techniques used by Kilcup et al. in zero temper- 
ature QCD fjl]. In addition to this and in order 
to help in the interpretation of our results we also 
computed the low lying spectrum of the Dirac op- 
erator ( in practice the lowest 8 eigenvalues and 
associated eigenvectors ) on each of our configura- 
tions. This was done using the conjugate gradient 
algorithm investigated by Kalkreuter and Simma 
B. As we will show later, some of the correlators 
are entirely saturated by the few lowest fermionic 
modes, allowing further insight into the dynamics 
of finite temperature QCD. 

All our computations were done on the CRAY 
C-90 at NERSC on lattices of size 16 3 x 8. We 
have gathered results both for quenched and 2 
flavor QCD. A single value of the quark mass has 
been studied so far: ma = 0.02 for the quenched 
case and ma — 0.00625 in 2 flavor QCD. In the 
later case, we use configurations generated by the 
HTMCGC collaboration g. The quenched com- 
putations are done at (3 = 5.8,5.9,6.0,6.1 and 
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6.2 where we have respectively analyzed 100, 100, 
100, 170 and 170 configurations. Similar num- 
bers for the 2 flavor case are (3 = 5.45 [160], 5.475 
[240], 5.5 [160], 5.525 [160] and 5.55 [80]. The 
position of the crossover is found to be slightly 
higher than (3 = 6.0 in the quenched case and 
slightly higher than (3 — 5.475 in two flavor QCD. 
In figure 1, we plot our results for the 7?, a and S 
screening lengths in 2 flavor QCD as a function 
of j3. The key feature of this plot is that the a 
becomes light close to the transition while the S 
remains heavy. This is in agreement with the ob- 
servation that the peak in the scalar susceptibility 
originates in the disconnected part of the corre- 
lator Above the phase transition, the Ua(1) 
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Figure 1. Screening masses in two flavor QCD 

symmetry breaking (m^ — ) appears to be of 
the same order as the explicit chiral symmetry 
breaking [m a — at our current value of the 
quark mass (ma = 0.00625). A careful extrapo- 
lation to the chiral limit will therefore be neces- 
sary in order to precisely sort out the situation 
in this region. This is currently being attempted 
by other groups although they only com- 

pute the integrated correlators (susceptibilities) 
instead of the full correlators as we do here. Using 
the data already available at this stage, however, 



it is tempting to speculate that the Uji(1) sym- 
metry is not restored within the critical region. 
Using the HTMCGC data for the <5 correlator at 
ma = 0.0125 § together with ours, it appears 
that mg extrapolates to a non-zero number at 
all values of (3. Assuming a second order tran- 
sition however, m a (and m^) will become very 
light close to the transition. Hence it appears that 
while SU(2) L x SU(2) R is being restored, U A (l) 
breaking is still present and only disappears at 
higher temperatures. 

As mentioned earlier we have also computed 
the lowest 8 eigenvalues and associated eigen- 
vectors on each of our configurations. Interest- 
ingly, we find that a spectral expansion of the 
quark propagator truncated to these low lying 
modes already captures a significant part of the 
physics. At or above the transition, for exam- 
ple, the pseudoscalar disconnected correlator is 
completely saturated by these low lying modes 
(similar computations below the phase transition 
are currently in progress). This suggests that a 
rather detailed description of the dynamics asso- 
ciated with the phase transition can be obtained. 
In the continuum, we expect to find two kinds 
of low lying modes: exact zero modes (which 
must appear in number consistent with the in- 
dex theorem: n+ — n_ = NfQ top ) and near zero 
modes (with no a priori connection to topology) . 
Modes from the first category satisfy r = 1 with 
r = | < n|75|n > |, whereas those from the second 
necessarily have r — 0. On the lattice however, 
these results are not reproduced exactly. Topo- 
logical zero modes are shifted away from zero and 
averages of the lattice T$ operator can differ sig- 
nificantly from the numbers given above §. By 
combining the two pieces of information (A and 
r), it is however often possible to decide on the 
topological or non-topological character of a given 
mode. This can be done for example by drawing 
a plot of | < n|75|n > | versus |A n |. Diagrams of 
this kind were considered by Hands and Teper in 
zero temperature QCD |ic|| . Here, we will focus 
our attention on the high temperature regime and 
leave the discussion of the transition region for a 
longer presentation . Our results for quenched 
QCD at j3 = 6.2 are gathered on figure 2. In order 
to check the consistency of our results with other 
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Figure 2. (r, A) diagram in quenched QCD at 
f3 = 6.2 



methods, we also cooled the configurations and 
assigned them an integer topological charge Q coo i- 
Results obtained on configurations with different 
Qcooi are represented by different symbols on fig- 
ure 2 (see insert). The results behave according 
to expectations; in particular, there is a strong 
correlation between large r and small eigenvalues 
and topological zero mode candidates can easily 
be identified. Note that they have a value of r 
around 0.20 only instead of 1 (This is to be ex- 
pected since T 5 being a 4-link operator picks up 
a large correction factor already in the mean-field 
approximation). By looking at these results con- 
figuration by configuration, one can check that 
the number of topological (zero) modes is exactly 
4:Qcooi (the 4 being associated with the 4 flavor 
of staggered fermions). The unique configuration 
with topological charge 2 in our sample, for ex- 
ample, gives us 4 low modes in figure 2 (the exis- 
tence of 4 other modes with opposite eigenvalue 
being guaranteed by the symmetries of the lat- 
tice Dirac operator). Finally, figure 2 also shows 
4 modes with low eigenvalues but small r. They 
all come from one given configuration with topo- 
logical charge 0. An analysis of the local values 
of < rt|r5(x)|n > reveals that this configuration 
contains an I — I pair and that the 4 fermionic 



modes are delocalized over the 2 topological ob- 
jects. In the continuum, this would give rise to a 
small but non-zero eigenvalue. On the lattice, the 
4 copies have different eigenvalues and different 
amount of derealization because of flavor sym- 
metry breaking. In conclusion of this analysis, 
we see that in quenched QCD at f3 = 6.2 the dis- 
connected pseudoscalar correlator is dominated 
by the contribution from topological zero modes 
(most of them with topological charge ±1). As 
the temperature is lowered, the number of small 
but non-exact zero modes increases and the re- 
gion close to the origin fills up in the (A, r) dia- 
gram. Finally, we should mention that the shift 
in the eigenvalue of topological modes (extending 
at least to the interval [-0.02,0.02] in figure 2) im- 
plies that we should be very careful in attempt- 
ing an extrapolation to the chiral limit. Using 
masses lower than 0.02 for example would lead 
to a significant underestimate of the quantities 
measured. It is also important to realize that 
the shift of topological modes is an ultraviolet 
effect (which could for example be corrected by 
the use of an improved or perfect action) and is 
therefore completely different from the depletion 
of eigenvalues of order 1/V commonly observed in 
measurements of < ipip > (which is a finite size 
effect). 
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